*Import industry share by location

use "$temp\Long_shares.dta" , clear

*Create industry information using data from ciccone and papaionau

gen capint=0.730 if ind39==22
replace capint=0.443 if ind39==24 
replace capint=2.103 if ind39==33
replace capint=2.316 if ind39==33
replace capint=0.481 if ind39==23
replace capint=1.746 if ind39==33
replace capint=1.744 if ind39==21
replace capint=3.194 if ind39==34
replace capint=1.366 if ind39==20
replace capint=2.215 if ind39==27
replace capint=1.807 if ind39==23
replace capint=0.785 if ind39==28
replace capint=2.265 if ind39==30
replace capint=0.789 if ind39==26
replace capint=1.173 if ind39==35
replace capint=2.385 if ind39==31
replace capint=1.632 if ind39==25
replace capint=1.199 if ind39==32
replace capint=1.320 if ind39==38
replace capint=1.017 if ind39==36
replace capint=0.800 if ind39==31
replace capint=0.924 if ind39==37

replace capint=0.654 if ind39==39
replace capint=1.416 if ind39==39
replace capint=1.954 if ind39==39
replace capint=0.878 if ind39==39


gen hcint=11.509 if ind39==22
replace hcint=11.244 if ind39==24
replace hcint=10.138 if ind39==33
replace hcint=11.547 if ind39==33
replace hcint=10.193 if ind39==23
replace hcint=13.204 if ind39==31
replace hcint=11.655 if ind39==32
replace hcint=11.9674 if ind39==21
replace hcint=11.435 if ind39==34
replace hcint=11.259 if ind39==20
replace hcint=11.693 if ind39==27
replace hcint=10.397 if ind39==23
replace hcint=12.792 if ind39==28
replace hcint=11.730 if ind39==30
replace hcint=10.760 if ind39==26
replace hcint=11.577 if ind39==35
replace hcint=12.704 if ind39==31
replace hcint=10.787 if ind39==25
replace hcint=11.921 if ind39==32
replace hcint=12.346 if ind39==38
replace hcint=11.484 if ind39==39
replace hcint=12.266 if ind39==36
replace hcint=13.031 if ind39==31
replace hcint=12.357 if ind39==37
replace hcint=12.518 if ind39==39
replace hcint=11.678 if ind39==39
replace hcint=11.354 if ind39==39


gen hcint1=hcint*share_ind
gen capint1=capint*share_ind

collapse hcint1 capint1, by(mun39)
rename hcint1 education_index1939
rename capint1 capital_intensity1939

save "$temp\education_capital_index_1939.dta", replace

